arXiv:1509.03908v2 [cond-mat.mtrl-sci] 12 Nov 2015 


Hybrid density functional study of optically active centers in GaN 

Khang Hoan£| 

Center for Computationally Assisted Science and Technology, 

North Dakota State University, Fargo, ND 58108, USA. 

(Dated: November 13, 2015) 

Understanding the luminescence of GaN doped with erbium (Er) requires a detailed knowledge 
of the interaction between the rare-earth dopant and the nitride host, including intrinsic defects 
and other impurities that may be present in the host material. We address this problem through a 
first-principles hybrid density functional study of the structure, energetics, and transition levels of 
the Er impurity and its complexes with N and Ga vacancies, substitutional C and O impurities, and 
H interstitials in wurtzite GaN. We find that, in the interior of the material, Erca is the dominant 
Er^"*" center with a formation energy of 1.55 eV; Eroa-VN possesses a deep donor level at 0.61 eV 
which can assist in the transfer of energy to the 4/-electron core. Multiple optically active Er®+ 
centers are possible in Er-doped GaN. 


1 Introduction Rare-earth dopants in semiconduc¬ 
tors have received great attention as they lead to sharp 
intra-/ shell optical transitions which are desirable for 
many optoelectronic applications. When optically ex¬ 
cited, Er-doped GaN, for example, will emit visible light 
in a few narrow bands in the green [l|. The Er cen¬ 
ter can be excited by a direct absorption of energy into 
the 4/-electron core (resonant excitation) or indirectly 
by energy transfer from the host (non-resonant excita¬ 
tion). The latter mechanism is believed to be mediated 
by Er-related defects in which the presence of defect lev¬ 
els in the host band gap, acting as carrier traps, may have 
a crucial role. It has been reported that multiple opti¬ 
cally active Er^"*' centers are present in the material; one 
center type may only be excited by resonant excitation 
while others are excited via non-resonant excitation 
[^. Song et al. observe different energy levels, at 0.188, 
0.300, 0.410, and 0.600 eV below the conduction-band 
minimum (CBM), presumably corresponding to different 
intrinsic and Er-related defects, in deep-level transient 
spectroscopy (DLTS) experiments on n-type GaN im¬ 
planted with Er 0. However, despite much effort, there 
is still no definitive assignment of these levels and the 
observed Er^+ centers to specific defect configurations. 

Computational studies of Er impurities in GaN have 
been carried out by several groups, e.g., using the local- 
density approximation (EDA) or self-interaction cor¬ 
rected LDA within density-functional theory (DFT) or 
the so-called LDA-l-1/ within a DFT-based tight-binding 
approach Filhol et ah, for example, report that 

substitutional Er at the Ga site (Erca) has a formation 
energy of 0.84 eV Q, while Sauna et al. find a much 
higher value (>6 eV) 0. These values are either much 
lower or higher than the experimental formation energy 
of 1.8±0.2 eV reported by Ug olini et al. for optically ac¬ 
tive Er^+ centers in GaN . We herein present a first- 
principles study of the Er impurity and its complexes 
with nitrogen and gallium vacancies, substitutional car¬ 
bon and oxygen impurities, and hydrogen interstitials, 
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using a hybrid functional approach. Although there are 
still limitations in the particular case of strongly corre¬ 
lated electrons, the method has been shown to be promis- 
ing in computational studies of defects in solids [Ti| . Car¬ 
bon, oxygen, and hydrogen are common unintentional 
dopants in GaN. In light of our results for defect struc¬ 
tures, formation energies, and transition levels, possible 
optically active Er^+ centers are identified. 

2 Methods Our calculations are based on DFT with 
the Heyd-Scuseria-Ernzerhof (HSE) functional [l^, as 
implemented in VASP [l3l |. The Hartree-Fock mixing 
parameter is set to 0.31, resulting in a band gap of 3.53 
eV for GaN, very close to the experimental value (^3.5 
eV). The defects are simulated using a 96-atom supercell, 
a plane-wave basis-set cutoff of 400 eV, and a 2x2x2 k- 
point mesh for the integrations over the Brillouin zone. In 
defect calculations, the lattice parameters are fixed to the 
calculated bulk values but all the internal coordinates are 
relaxed; spin polarization is included. The Ga 3d and Er 
4/ electrons are treated as core states within the projec¬ 
tor augmented wave method MM- Lyons et al. reports 
that the inclusion of the Ga semicore d electrons in the 
valence has negligible effects on the formation energies 
and transition levels of defects in GaN [T^. Our HSE 
calculations of ErGa 47 N 48 in which the Er 4/ electrons 
are explicitly considered show that there are no / states in 
the band gap or near the valence-band maximum (VBM) 
and CBM, indicating that the explicit inclusion of the / 
electrons in the valence would have negligible effects on 
the energetics of Er-related defects in GaN. Note that 
un-spin-polarized DFT-|-(7 results reported by Wang et 
al. [13 also show the VBM and CBM are predominantly 
bulk GaN states. 

The likelihood of a intrinsic defect, impurity, or defect 
complex X in charge state q being incorporated into a 
crystal is characterized by its formation energy, defined 
as 

E;-^(X9) = E;tot(X«) - Atot(bulk) (1) 

i 

+ q{Ey A- He) -b A'^, 

where EtotO^'^) and £'tot(bulk) are the total energies of 
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FIG. 1. Formation energies of Vn, Vgh, Erca, and Eroa- 
Vn and Erca-Fba complexes in GaN, plotted as a function of 
Fermi level, under Ga-rich and N-rich conditions. Fermi levels 
above the CBM (at 3.53 eV) are indicated by the shaded area. 


the defect and bulk supercells; rii is the number of atoms 
of species i that have been added (rn > 0) or removed 
{rii < 0) to form the defect; fii is the atomic chemical 
potential, representing the energy of the reservoir with 
which atoms are being exchanged. The chemical poten¬ 
tials of Ga, C, N, and H are referenced to the total energy 
per atom of bulk Ga, bulk G (diamond), N 2 and H 2 at 
0 K, while those of Er and O are assumed to be limited 
by the formation of ErN and Ga 203 . The formation en¬ 
thalpies {AH, at 0 K) are calculated to be —1.26, —4.00, 
and —10.07 eV for GaN, ErN, and Ga 203 , in agreement 
with the experimental standard formation enthalpies of 
-1.63±0.17 [13, -3.71±0.23 [H, and -11.20 eV [13. 
The chemical potentials of Ga and N can vary over a 
range determined by the formation enthalpy of GaN such 
that ^Ga + Mn = Ai7(GaN); = 0 (/tn = 0) corre¬ 
sponds to Ga-rich (N-rich) condition, /ig is the electronic 
chemical potential, i.e., the Fermi level, representing the 
energy of the electron reservoir, referenced to the VBM 
in the bulk (Ev). A'^ is the correction term to align the 
electrostatic potentials of the bulk and defect supercells 
and to account for finite-size effects on the total energies 
of charged defects [2ll |. 

From defect formation energies, one can calculate ther¬ 
modynamic transition levels (not to be confused with op¬ 
tical transition levels) between charge states q and q' of 
a defect, e{q/q'), defined as the Fermi-level position at 
which the formation energy of the defect in charge state 
q is equal to that in charge state q' [ll|. e{q/q') is in¬ 
dependent of the atomic chemical potentials. This level 
would be observed in, e.g., DLTS, where the defect in the 
final charge state fully relaxes to its equilibrium configu¬ 
ration after the transition. 

3 Results and discussion In bulk GaN, each Ga 
atom is coordinated with four N atoms: one along the 


c-axis with the Ga—N bond length calculated to be of 
1.958 A and three in the basal plane with the Ga—N 
bond length of 1.952 A. Figure [1] shows the calculated 
formation energies of the nitrogen vacancy (I4 j), gallium 
vacancy (Vca), Eroa, and complexes between Ei'ca and 
the vacancies. Vn is found to be stable in the 3-I-, -I-, and 
neutral charge states; the axial and basal positions of Vn 
are degenerate in energy. This vacancy introduces two 
transition levels in the band gap: one at 0.53 eV above 
the VBM corresponding to the transition between the 3-1- 
and -I- charge states and the other at 0.27 eV below the 
CBM corresponding to the transition between the -I- and 
neutral states, in good agreement with those reported by 
Yan et al. [^. Vca is stable in all charge states from 
3— to +. The (-I-/0) transition level occurs at 1.03 eV, 
followed by the (0/—) level at 1.67 eV, the (—/2—) level 
at 2.24 eV, and the (2-/3-) level at 2.82 eV above the 
VBM, in agreement with Lyons et al. [^ . 

We find that Erca is stable only in the neutral state, 
c/. Fig.[U This is also observed in our calculations where 
Er 4/ electrons are included in the valence. The Er^+ 
ion does not only introduce an outward relaxation of the 
neighboring N ions, but also distort the local lattice envi¬ 
ronment, making the trigonal distortion much more pro¬ 
nounced with the Er—N bond length along the c-axis is 
~0.03 A longer than the other bonds; the Er—N bond 
lengths are calculated to be 2.141, 2.143, and 2.148 A in 
the basal plane and 2.169 A along the c-axis. Er^+ also 
goes off-center by moving toward the basal N plane by 
0.08 A. These results are in contrast to previous work 
which argues that Er^+ remains “exactly on-center” and 
even reduces the difference in the axial and basal bond 
lengths The formation energy of Erg^^ is cal¬ 

culated to be 1.55 eV and remains unchanged in going 
from Ga-rich to N-rich conditions, as it depends on the 
atomic chemical potentials only through the difference in 
the formation enthalpies of GaN and ErN. This value cor¬ 
responds to, e.g., a concentration of 5x10^® cm“^ under 
thermal equilibrium at 1040°C. In the absence of defects 
and other impurities in its close proximity, Erca would 
thus be electrically inert. Still, intra-4/ shell optical tran¬ 
sitions can occur in these isolated Er^+ centers as the lo¬ 
cal lattice distortion mentioned above breaks the inver¬ 
sion symmetry, thus relaxing the parity selection rules. 

Regarding Erca-Wj, there are two configurations cor¬ 
responding to different positions of Vn in the complex: 
“axial” or “basal” as shown in Fig. [2] They both have 
large local lattice distortions; Er^+ moves off-center to¬ 
ward the vacancy by 0.22 A (0.20 A) and the distance 
between the two components in (Erca-VN)® is 1.73 A 
(1.77 A) in the basal (axial) configurations. The basal 
configuration has a lower formation energy, as seen in 
Fig. [T] In the following, we refer to a defect complex 
only in its lowest-energy configuration. Note that the lo¬ 
cal distortion is also different for different charge states. 
In (Ei'Ga-lAi)^’'", for example, the distance between the 
two constituents is 1.66 A; Er^+ thus moves toward the 
vacancy by 0.29 A. The binding energy of (ErGa-VN)'^, 
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FIG. 2. Structure of (ErGa-Vk)'^: (a) axial and (b) basal 
configurations. Large circles are Er, medium Ga, and small 

N. 


with respect to its isolated constituents Eiq^ and is 
0.94, 0.60, 0.74, 0.56, or 1.28 eV for q = 0, +, 2+, 3+, 
or —, respectively. The binding energy is thus relatively 
high (0.94 eV) in n-type GaN where Erca-^N is most 
stable in the neutral state, cf. Fig. [U suggesting that a 
sig nificant portion of Vn can exist in form of Erca-EN 
[ 3 ]. This defect complex has several deep levels in the 
band gap, corresponding to the (3+/2+) level at 0.38 eV 
and the (2+/+) level at 0.64 eV above the VBM, and the 
(+/0) level at 0.61 eV below the CBM. The (0/—) level is 
only 0.02 eV above the CBM, cf. Fig. [TJ The deep donor 
level at 0.61 eV below the CBM can act as a carrier trap; 
an electron captured here can subsequently recombine 
nonradiatively with a free hole from the valence band or 
a hole at some acceptor level and transfer energy to the 
4/ shell of the Er^+ ion. We note that the concentra¬ 
tion of Erca-f^ should be much smaller than that of the 
isolated Erca! e.g., assuming n-type GaN and extreme 
Ga-rich condition, the formation energy of the complex 
is 3.82 eV, corresponding to an equilibrium concentra¬ 
tion of only 3x10® cm“® at 1040°C. However, the cross 
section for electron capture associated with the complex 
is expected to be larger, as indicated in experiments Q. 
Note that the concentration of certain defects could be 
much higher near the surface or if the material is pre¬ 
pared via a non-equilibrium process. 

For (ErGa-Vca)'^, we find a binding energy of 0.44, 0.26, 
0.24, 0.15, and 0.02 eV with respect to ErQ,^ and for 
q = 3—, 2—, —, 0, or +, respectively. These values are 
thus much smaller than those in the case of (ErGa-f^)"^, 
mainly due to the larger distance between ErGa and Vca- 
For example, the distance between the two components 
in (ErGa-Vca)^”, the most stable charge state in n-type 
GaN, is 2.75 A. ErGa-Vca induces several deep levels in 
the band gap: the (-b/O), (0/—), (—/2—), and (2-/3—) 
levels occur at 0.91, 1.58, 2.22, and 2.65 eV above the 
VBM, respectively. Some of these levels may act as car¬ 
rier traps and thus assist in the transfer of energy to 
the Er 4/ shell. However, given the complex’s relatively 
low binding energy and high formation energy (e.g., com¬ 
pared to ErGa-lAr, except at the Fermi levels very close to 
the CBM and under extreme N-rich condition, cf. Fig.[T]), 
the role of ErGa-Vca in the Er-related luminescence is still 
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FIG. 3. Formation energies of On, Gn, and Hi, and their 
complexes with Erca in GaN under Ga-rich and N-rich con¬ 
ditions. Fermi levels above the CBM are indicated by the 
shaded area. 


not clear at this point. 

Figure [3] shows the calculated formation energies of the 
substitutional oxygen (On) and carbon (Cn) at the N site 
and hydrogen interstitial (Hi), and their complexes with 
ErGa- On is found to be a shallow donor. Cn induces two 
levels in the band gap: the (-F/O) level at 0.31 eV and the 
(—/O) level at 1.02 eV above the VBM. Hi is amphoteric, 
being positively charged in p-type and negatively charged 
in n-type GaN; the (-1-/—) level occurs at 3.03 eV above 
the VBM. Our results for On, Cn, and Hi are thus in 
good agreement with those reported previously [H, [H, 

13 . 

The energy landscape of the complexes between ErGa 
and On, Cn, and Hi resembles that of the isolated non-Er 
impurities, cf. Fig.[3l Like On, ErGa-ON is also a shallow 
donor; the complex in its -F charge state has a binding 
energy of 0.60 eV with respect to Ecq^^ and OJ^. Since 
there are no in-gap levels associated with ErGa-ON, this 
complex is unlikely to be effective in assisting the energy 
transfer to the Er 4/ shell. ErGa-CN induces two levels in 
the band gap: the (-F/0) level at 0.11 eV and the (0/—) 
level at 2.05 eV. However, the most stable charge state of 
the complex in n-type GaN, (ErGa-CN)“, has a negative 
binding energy of —0.63 eV and is thus unstable with 
respect to Er^,,^ and Cj;^. The binding energy of other 
charge states is also quite low: 0.40 or 0.20 eV for the 
neutral or -F state. We also consider the ErGa-ON-CN 
complex where O is at one of the basal N positions and C 
is at the axial position; the binding energy is 0.49 eV in its 
neutral state or 0.94 eV in the -F state. The (-F/0) level 
induced by this complex is, however, quite low at 1.46 eV 
above the VBM, cf. Fig.[3l thus excluding its possible role 
in high-energy luminescent transitions. Finally, ErGa-H^ 
is found to have a binding energy of only 0.04 or 0.25 eV 
in the -F or — charge state, suggesting that it is unlikely 
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to be stable as a defect complex. 

In n-type GaN, is spin-polarized with spin S = 1/2 
and (ErGa-14^)° has a total spin of 2 as seen in our calcu¬ 
lations where the Er 4/ electrons are explicitly included; 
the component of (Erca-I^i)” has S = 1. Other 
defects relevant to n-type GaN are not spin-polarized, 
except the Er^^^^ component which has spin S = 3/2. Ex¬ 
perimentally, Ugolini et al. obtain a formation energy of 
1.8±0.2 eV for Er^+ centers in Er-doped GaN prepared 
by MOCVD. These authors presume the centers to be 
ErGa-f4j 0- We find that, given the much higher forma¬ 
tion energies of both En and the complex in n-type GaN, 
cf. Fig. [U those Er^+ centers should instead be identi¬ 
fied with isolated ErGa which has a calculated formation 
energy of 1.55 eV. Regarding En, the (-I-/0) transition 
level can be identified with the DLTS level at 0.270 eV 
in undoped GaN or 0.300 eV in Er-implanted GaN ob¬ 
served in experiments Among the remaining DLTS 
levels reported by Song et ah, the level at 0.600 eV can 
be identified with the (-I-/0) level of ErGa-f4j; other levels 
at 0.188 eV (often mentioned in the literature as being 
related to ErGa-f4j) and 0.410 eV are likely to involve 
some other defect structures including intrinsic and/or 


Er-related defects near the surface. Further studies are 
needed to clarify the defect physics of Er-doped GaN. 

4 Summary and conclusions In summary, we find 
that isolated ErGa is the dominant Er^+ center, at least 
in n-type GaN where it has the lowest formation energy 
among all Er-related defects. This impurity possesses 
a local lattice distortion that can help relax the selec¬ 
tion rules, thus allowing for intra-4/ transitions to occur 
in the Er^+ ion. ErGa-^N is identified as an efficient 
defect-related Er center for band-to-band (non-resonant) 
excitation. This complex has the highest binding energy, 
shortest distance between the defect and Er, and a deep 
trap at 0.61 eV below the CBM. ErGa-Vca, in principle, 
can also act as an assistant for the transfer of energy to 
the Er 4/ shell. Multiple luminescent centers with, e.g., 
different Stark splittings, luminescence decay dynamics, 
and excitation cross sections are thus possible in real Er- 
doped GaN samples. 
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